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We propose to study hadron-hadron scattering using lattice QCD in an asymmetric box which 
allows one to access more non-degenerate low-momentum modes for a given volume. The con- 
ventional Liischer's formula applicable in a symmetric box is modified accordingly. To illustrate 
the feasibility of this approach, pion-pion elastic scattering phase shifts in the 7 = 2,7 = channel 
are calculated within quenched approximation using improved gauge and Wilson fermion actions 
on anisotropic lattices in an asymmetric box. After the chiral and continuum extrapolation, we 
find that our quenched results for the scattering phase shifts in this channel are consistent with the 
experimental data when the three-momentum of the pion is below 300MeV. Agreement is also 
found when compared with previous theoretical results from lattice and other means. Moreover, 
with the usage of asymmetric volume, we are able to compute the scattering phases in the low- 
momentum range (pion three momentum less than about 350MeV in the center of mass frame) 
for over a dozen values of the pion three-momenta, much more than using the conventional sym- 
metric box with comparable volume. 
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1. Introduction 

Hadron-hadron scattering experiments play a crucial role in understanding interaction among 
hadrons and QCD. Due to its non-perturbative feature, lattice QCD (LQCD) is the only systematic, 
non-perturbative method of QCD which in principle can be applied to calculate low-energy hadron- 
hadron scattering from first principles. This type of lattice calculations is realized by a finite size 
method proposed by Liischer [jTJ]. Using this technique, the scattering length and the scattering 
phase shifts for pion-pion scattering in the I = 2, J = channel have been studied P, |L 0, |L ||, [7], 
H|, |9, [lO, O, 12] in both quenched and unquenched lattice QCD. 



However, non-degenerate low-momentum modes accessible for lattice simulation are usually 
limited due to computational cost. With cubic boxes, more non-degenerate low-momentum modes 
requires larger physical volumes. To circumvent this problem, we propose to use asymmetric boxes 
to study hadron-hadron scattering. The idea is tested in a quenched study of pion-pion scattering 



in the / = 2, J = channel [24] using improved anisotropic lattices. 



2. Liischer 's formulae extended to an asymmetric box 

Consider a cubic box with size LxLxL and periodic boundary condition. Three momentum 
of a single pion is quantized as: k = (2n/L)n with n = (ni,ri2,nj) € Z 3 , where Z represents the 

-2 

set of all integers. For the two-pion system taking the center of mass reference frame, we define k 
of the pion pair in a box as: 



E KK = 2^ml + \L (2.1) 

where E n% is the exact energy of a two-pion system with the two pions having three momentum k 
and — k respectively in the center of mass frame. We further define q via: 

k 2 = ^ 9 2 (2-2) 
With these definitions, Luscher's formula for the s-wave scattering phase can be represented as: 

cot5o® = ^4^ 1 ' ( 2 - 3 ) 



where the so-called zeta function 3fi m is given by: 



^sd) = Y.T=F^L- ( 2 - 4 ) 



In this definition, %,(r) = r'F/ m (£2 r ), with Yi m (Q. r ) being the usual spherical harmonics. 

We need to extend the above formulae to an asymmetric box of physical size: L x 772L x T73L. 
The momentum in the box is quantized as: k = (2n/L)n with n = (rai,^/^;^/^) an d n = 
(«i,«2i"3) S Z 3 . Quantities k and q are still defined according to Eq. (2.1) and Eq. ( |2.2| ). The 
formula for the scattering phase shifts is now modified to: [ |13i Q 



t * n\ 1 \ — ^Oo(l, <7 .172,173) n n 
cotdo(l) = m m (q) = T75 (2-5) 
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with the modified zeta function 3?i m defined as: 



%m (S, q 2 ; TJ2 , TJ 3 ) = £ 



n - Q 

In the large volume limit, the formula for scattering length is: 



2Vv 



2m % 



l+Ci(Tj2,TJ 3 ) 



+ C 2 



<*-*>(r) ! 



+ 0(ZT 6 ). 



(2.6) 



(2.7) 



where the coefficients for the case T] 2 = 1 and T] 3 = 2, which is what we simulated in this work, are 
found to be: ci(l,2) = -1.805872, c 2 (l,2) = 1.664979. 

3. Simulation details and results 

To test our idea of using the asymmetric box on hadron-hadron scattering, we perform a 
quenched study on the pion-pion scattering phase shift in the / = 2, J = channel. In this sec- 



tion, we will briefly introduce our numerical results. More details can be found in Ref . [ 24 ] 



p 


Lattice 


a s (GeV- 1 ) 


Confs. 


2.080 


8 2 x 16 x 40 


1.5677 


464 


2.215 


9 2 x 18 x 48 


1.3926 


425 


2.492 


12 2 x 24 x 64 


1.0459 


105 



The gauge action used in this study is the tadpole improved gluonic action on anisotropic lat- 
tices [|l(], 17] and the fermion action used in this calculation is the tadpole improved clover Wilson 
action on anisotropic lattice [18, 19, 2(], 21, |2|]. The anisotropy parameter E, = a s /a t is fixed to 5 
for all our lattices. Some of the simulation parameters are provided in the table above. 

To obtain the two-pion energy levels, we follow the standard procedure of correlation matrix 
measurements. For the single pion operators, we use: 



x n 



-d(x,t)r 5 u(x,t), 



(3.1) 



The Fourier transformed fields are defined accordingly. The s-wave two-pion operators in the 1 = 2 
channel are defined as: 

^) = ^ (K Jt)^._ K) (t + l). (3.2) 



where the summation runs over all symmetry operations R of the symmetry group for the asym- 
metric box. In the simulation, we measure the correlation matrix among different non-degenerate 



two-pion modes, using the two-pion operators defined in Eq. (3.2): 

%*(t) = (0i(t)0 n (t s )) . 

We then constructed a new correlation matrix: 

a{t,t Q ) = <g(t Q )-itf(t)tf(t Q )--i 

where to is some suitable reference time. The eigenvalue of this new matrix £2 is: 



(3.3) 



(3.4) 



(3.5) 
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Figure 1: Effective mass plateaus of two-pion system at j3 = 2.492 after the diagonalization procedure. The 
horizontal line segments in the figure represent the fitting ranges of the plateaus. Results for other values of 
j3 are similar. Only two momentum modes out of 6 are shown. 

from which the two-pion energy eigenvalues Ei are obtained. In Fig. [I], we have shown the effective 
mass plateaus for the two-pion energies at j3 = 2.492. The horizontal line segments designate the 
fitted values for the plateaus. Due to the usage of anisotropic lattices, good accuracy is achieved 
for two-pion energy levels. 

In our quenched study, the pion mass values are beyond the applicable range of chiral perturba- 
tion theory. We therefore adopted polynomial parametrization for the data points for the scattering 
lengths and the phase shifts. For the scattering length, we parameterize our data for ao / m n using 
either a linear function in nv\: 

^=A + A l m 2 7t ;, (3.6) 
or, for the case of j8 = 2.492 where a curvature is seen, a quadratic form in m„. After the chiral 




Figure 2: Left panel: chiral extrapolation of the quantity Right panel: the corresponding continuum 
limit extrapolation of the quantity The extrapolated result is indicated by a solid circle near a s = 0. 

extrapolation, the results for the scattering length at different lattice spacing values are extrapo- 
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lated towards the continuum with a linear function in a 2 . The chiral extrapolation and the con- 
tinuum extrapolation for the scattering length are shown in Fig. |2[ As a result, we finally obtain: 
aom n = 0.0359(59) which is compatible with results obtained in other lattice simulations, chiral 
perturbation theory and the experiment. 

Similarly, we parameterize the phase shift 8 itself by a polynomial in both m 2 and k 2 as: 



S(ml,k 2 ) = D m + D w m z n + D 2 ml + Do X k 1 +D n m 1 K V+D ()2 k 



2 7,2 



(3.7) 



The coefficient Dqo vanishes due to chiral symmetry which is checked numerically as well. A linear 
extrapolation in a 2 s is then performed for various coefficients Djj appearing in the above equation. 
Some of the extrapolations are shown in Fig. ||[ 
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Figure 3: The continuum limit extrapolation of the parameter Dy. The solid circles near a s — represent 
the corresponding continuum limit values. 

The results for the phase shift after the chiral and continuum extrapolations are summarized 
in Fig. ^ together with the experimental results from CERN-Munich group Jl5|]. First, it is readily 
seen that, due to the usage of the asymmetric box, much more low-momentum modes are accessible 
with a relatively small volume, allowing much better comparison with the experimental data in the 
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Figure 4: Comparison of our lattice results for the scattering phase shifts with the experimental data from 
CERN-Munich [|l5|]. Results are consistent with each other for k 2 below 0.1 GeV 2 which roughly corre- 
sponds to the center of mass energy of about 0.6GeV. 



low-momentum range. Second, it is also found that our final results agree with the experimental 
results within errors for k below 0.1 GeV which is about y s = 0.6 GeV. At higher energies, our 
results deviate from the experimental results. This deviation might be caused by the systematic 
uncertainties in our calculation, e.g. quenching and chiral extrapolations. Numerically speaking, 
it is largely due to poor determination of the coefficient D02. However, we would like to point out 
that, it is clear from our quenched calculation that, the asymmetric volume technique advocated 
here would also be useful for unquenched studies once the unquenched configurations become 
available. 



4. Conclusions 

In this paper, we propose to study hadron-hadron scattering processes on lattices with asym- 
metric volume. This setup has the advantage that it provides much more non-degenerate low- 
momentum modes with a relatively small volume, allowing more detailed comparison both with 
the experiments and with other theoretical results. The feasibility of this idea is illustrated in a 
quenched calculation of pion-pion scattering length and scattering phases in the I = 2, J = chan- 
nel using clover improved lattice actions on anisotropic lattices. Our quenched results indicate that 
the usage of asymmetric volumes indeed allow us to access much more low-momentum modes than 
in the case of cubic volume of similar size. For k 2 in the range of 0.02GeV 2 to 0.12GeV 2 , we have 
over a dozen of data points for the phase shift, much more than that in the cubic case with similar 
volume. It is also noted that, in the low-momentum region, after the chiral and continuum extrap- 
olations, our results for the scattering length and the scattering phase shifts are in good agreement 
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with the experimental data and are consistent with results obtained using other theoretical means. 
Although our calculation is now performed in the quenched approximation, similar calculations are 
also possible in the unquenched case once the gauge field configurations become available. 
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